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^ , Abstract 



We study the origin and bifurcations of typical classes of unstable periodic orbits in a 
Q , jet flow that was introduced before as a kinematic model of chaotic advection, transport 

I and mixing of passive scalars in meandering oceanic and atmospheric currents. A method 

' to detect and locate the unstable periodic orbits and classify them by the origin and 

bifurcations is developed. We consider in detail period-1 and period-4 orbits playing an 
important role in chaotic advection. We introduce five classes of period-4 orbits: western 
and eastern ballistic ones, whose origin is associated with ballistic resonances of the fourth 
order, rotational ones, associated with rotational resonances of the second and fourth or- 
ders, and rotational-ballistic ones associated with a rotational-ballistic resonance. It is a 
new kind of nonlinear resonances that may occur in chaotic flow with jets and/or circula- 
tion cells. Varying the perturbation amplitude, we track out the origin and bifurcations 
of the orbits for each class. 

> 

o ■ 1 Introduction 

fSl ■ It is well known that a dynamical system is chaotic if it displays sensitivity to initial conditions, 

has a dense orbit and a dense set of periodic orbits (see, for example, [1]). Periodic orbits play 
Q ! an important role in organizing dynamical chaos both in Hamiltonian and dissipative systems. 

Stable periodic orbits (SPO) organize a regular motion inside islands of stability in the phase 
space. Unstable periodic orbits (UPO) form a skeleton around which chaotic dynamics is 
organized. The motion nearby an UPO is governed by its stable and unstable manifolds. 
Owing to the density property, the UPOs influence even the asymptotic dynamics. Order and 
disorder in a chaotic regime are produced eventually by an interplay between sensitivity to 
initial conditions and regularity of the periodic motion. 

In the present paper we study origin and bifurcations of typical classes of the UPOs in a two- 
dimensional incompressible flow that has been introduced and analyzed in Refs. [21 [Sj HI O [6] as 
a toy kinematic model of transport and mixing of passive particles in meandering jet currents 
in the ocean, like the Gulf Stream, Kuroshio and other main oceanic currents, and in the 
atmospheric currents like zonal jets with propagating Rossby waves fl2 [ [H I [16]. The equations 
of motion of passive particles advected by any incompressible planar flow are known to have a 
Hamiltonian form 

^ = ^x,y,t) = -— , ^t=vi^,y^t) = ^, (1) 

where the streamfunction \E' plays the role of a Hamiltonian, and the particle's coordinates 
X and y are canonically conjugated variables. The phase space of Eqs. [1] is a physical space 
for advected particles. If the velocity field, u = u{x,y) and v = v{x,y), is stationary, then 
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Figure 1: The phase portrait of the model flow ([3]) in the comoving frame of reference. The 
first frame with streamlines in the circulation (C), jet (J) and peripheral currents (P) zones is 
shown. The parameters are: A = 0.785, C = 0.1168 and L = 0.628. 

fluid particles move along streamlines, and the motion is completely regular with any Eulerian 
stationary field whatever its complexity. A time-periodic velocity field, u{x, y, t) = u{x, y,t + T) 
and v{x, y, t) = v{x, y, t+T), can produce chaotic particle's trajectories, the phenomenon known 
as "chaotic advection" [TllH]. 

Chaotic advection of water (air) masses along with their properties in geophysical jets is 
a topic of great interest in the last decade (for recent reviews on chaotic advection, transport 
and mixing in the ocean and atmosphere see [9] and [10], respectively). Among the variety 
of kinematic and dynamic models of shear flows, one of the simplest ones is a Bickley jet 
with the velocity profile ~ sech^?/ and a running wave imposed. The phase portrait of such 
a flow in the frame, moving with the phase velocity of the running wave, is shown in Fig. [H 
The flow consists of three distinct regions, the eastward jet (J), the circulations (C) and the 
westward peripheral currents (P) to the north and south from the jet, separated from each 
other by the northern and southern oo-like separatrices. A simple periodic modulation of the 
wave's amplitude breaks up these separatrices, produces stochastic layers in place of them, and 
chaotic mixing and transport of passive particles may occur. 

In the recent papers [H [6] we have studied statistical properties of chaotic mixing and 
transport in such a time-periodic meandering Bickley-jet current and explained some of them by 
the presence of dynamical traps in the phase space, singular zones in the stochastic layers where 
particles may spend arbitrary long but finite time [H]. We identified rotational-islands traps 
around the boundaries of rotational islands, ballistic-islands traps, around the boundaries of 
ballistic islands, and saddle traps associated with stable manifolds of periodic saddle trajectories 

A further insight into chaotic advection in the chosen model flow and in other jet flows is 
required in order to find a connection between dynamical and topological properties of particle's 
trajectories forming a complex picture of chaotic mixing. The aim of this paper is to study in 
detail origin, properties and bifurcations of typical UPOs in the flow considered in Refs. [HE]- In 
Sec. 2 we introduce the model streamfunction and advection equations and present a numerical 
method for locating the UPOs of different periods in chaotic dynamical systems. The method is 
based on computing a distance d between the positions, x(to), yito), of a chosen particle at the 



moments of time to aiid to + (where To is a period of the perturbation and m = 1, 2, . . . ), 
finding local minima of the distance function with given values of m and analyzing them to 
locate the period-m UPOs for which d{x{to), ?/(to)) = 0. In Sec. 3 we analyze a saddle orbit (SO) 
by linearizing the advection equations and study its metamorphoses varying the perturbation 
amplitude e. In Sec. 4 we apply the method to study the origin, properties and bifurcations 
of the period-4 UPOs. We chose namely that period because for sufficiently large values of 
e there are no visible period-4 resonances in the phase space. It is easy to locate the UPOs 
and SPOs of visible resonances on Poincare sections, but it is not a trivial job to do that with 
broken resonances. All the period-4 UPOs are classified in 5 distinct groups which differ by 
the type of motion of the corresponding particles and their origin and bifurcations. Varying 
the perturbation amplitude e, we search for those resonances in each class which generate the 
corresponding period-4 UPOs and compute bifurcation diagrams for them. 



2 Return maps for the model flow 

Following to Refs. |H |6], we consider the following streamf unction in the laboratory frame of 
reference: 

, , ,, ^, , / y' — a cos k(x' — ct') 

\X^l + k^a^ smH{x' -ct' 

where the hyperbolic tangent produces the Bickley-jet profile, the square root provides a con- 
stant width of the jet A, and a, k and c are amplitude, wave number and phase velocity of the 
running wave, respectively. The normalized streamfunction in the frame moving with c is 

^ / y — A cos X \ ^ 
^ = -tanh ^ ^ = + Cy, (3) 



LVI + sin^ 



X 



where x = k{x' — ct') and y = ky' are new scaled coordinates. The normalized jet's width L = 
Xk, wave's amplitude A = ak and phase velocity C = c/'^^k are the control parameters. The 
advection equations ([T]) with the streamfunction ([3]) have the following form in the comoving 
frame: 

1 . Asm. x{l + A^ — Ay cos x) 



h 



L Vl + A^ sin' xcosh^ d ' L (l + ^2 sin'a;)^^' cosh' 0' 

y — A cos X 



(4) 



L Vl + A^ sin' X 

where dot denotes differentiation with respect to the scaled time t = '^'^k^t' . 

The flow with the streamfunction ([2]) is steady in the comoving frame and its phase portrait 
is shown in Fig. [TJ There are southern and northern sets of elliptic fixed points: Xe = 27m, 
yi^^ = —L Arcosh -^XjLC + A and Xe"^ = (2n + l)7r, ?/e"^ = L Arcosh ^Jl/ LC — A, respectively, 
and the southern and northern sets of hyperbolic (saddle) fixed points: xi'^'* = {2n + l)7r, 
yi^"* = —L Arcosh ^\jLC — A and x^f^ = 27m, yi^'^ = L Arcosh ^JXjLC + A, respectively, 
where n = 0, ±1, .... 

A perturbation is provided by a periodic modulation of the wave's amplitude 

A{t) = Aq^ ecos{ujt^ Lp). (5) 

The equations of motion (jl]) are symmetric under the following transformations: t —>■ t, x ^ 
71 + X, y ^ —y and t —t, x —x, y ^ y. Due to these symmetries, the motion can 
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be considered in the northern chain of the circulation cells on the cylinder with < a; < 27r. 
The part of the phase space with 27m < x < 2TT{n + 1), n = 0, ±1, . . . , is called a frame. The 
first frame is shown in Fig. [H The values of the following control parameters are fixed in our 
simulation: L = 0.628, = 0.785, C = 0.1168, Tq = 27r/cu = 24.7752, ip = 7r/2. The only 
varying parameter is the perturbation amplitude e. 

In fluid mechanics an infinite number of initial conditions comes in play simultaneously and 
a number of fluid elements, launched in different places, may follow the same orbit on the flow 
plane. Essentially, that a number of particles with different initial positions may move along 
the same orbit. More precisely, an orbit is a set of points Xi,yi [i = 1,2,...) on the phase 
plane (on the flow plane) with the following two properties: (i) there exists for Vz, j an integer 
k (positive or negative) such that 



where U is an evolution operator; and (ii) there exists for Vi, /c an integer j such that 



A period-m orbit is a finite set of points on the phase plane with the properties ([6]) and (jT]) 
consisting of m elements. Thus, any period-m orbit contains m points whose trajectories belong 
to this orbit. 

To locate the UPOs in the phase space we fix values of the control parameters and compute 
with a large number of particles the Euclidean distance d"^ = [x{to + mTo) — x(to)]^ + [yitQ + 
mTo) — ?/(to)]^ between particle's position at an initial moment of time to and at the moments of 
time T = mTo, where m = 1, 2, . . . . The data are plotted as a period-m return map (RM) that 
shows by color the values of d for particles with initial positions [xito], yito)]. At the first stage, 
we select a large number of points where the function d{x{to),y{tQ)) may have local minima. 
Then we apply the method of a deformed simplex to localize the minima in neighbourhoods 
of those points. There are such minima among them for which d = with a given value of 
m. The procedure allows to detect both UPOs and SPOs not only in periodically perturbed 
Hamiltonian systems but in any chaotic system. 

Return maps with m = 1, 2, 3, 4 and 12 have been computed. In the next section we consider 
the SO which is a period-1 UPO. The main efforts are devoted to analysis of period-4 UPOs 
because it is not a trivial task to detect the UPOs when the corresponding resonances cannot be 
identified on Poincare sections. It is the case with m = 4. The period-4 RM is shown in Fig. [2^ 
with the cross marking location of the period-1 saddle trajectory and black dots marking initial 
positions of twenty six trajectories of the period-4 UPOs with a relative accuracy 10~^^ -h 10~^^. 
For comparison, we demonstrate in Fig. [2)3 locations of those dots on the Poincare section. It 
is evident that they cannot be prescribed to any structures in the phase space and cannot be 
identified by inspection of the Poincare section. 

3 Saddle orbit 

The saddle points of the unperturbed equations of motion (jl]), (xi""*, yi"^) and {xi^\yi^^), become 
period-1 saddle orbits (SO) under the periodic perturbation ([5]). To analyze such a SO we 
linearize the perturbed equations 
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x = X{x,y), y = Y{x,y) 
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Figure 2: (a) Period-4 return map representing the distance d between particle's position x{to) 
and y(to) at to and its position at to + 4To. Cross marks location of the saddle trajectory and 
dots mark initial positions of trajectories of the period-4 UPOs. (b) Poincare section of the 
northern separatrix layer with positions of those dots. 
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in a neighborhood of the saddle point 



Xs = 0, ys = L Arcosh ^ + ^o, (9) 

where X{x, y) and Y{x, y) are the right-hand sides of the corresponding equations in the set (jlj) 
with A being the time-periodic amphtude ([5]). After hnearizing, we get the following equations 
for small deviations rj = x — Xg and ^ = y — yg-. 

1 ^ 2tanhe ^ • A{ALQ-l) 

V = 9 C ^ = — ^ — -ri, (10) 

' Lcosh^e L2cosh2e Lcosh^e ' ^ ' 

where 



^ . , / 1 £:cos$ 
e = Arcosh W— — , $=cjt + v?. 

For small values of the perturbation amplitude e, it is possible to simplify the set ffTOl) 
expanding the right-hand sides in a series in powers of e cos $ and neglecting terms above the 
first order. In terms of the variables 

Xs = limX{xs,ys), Xys = hmXy{xs,ys), = limK,(x5, 

£^0 e^O e^O 



6^0 d{e COS $) ' £-0 d{e cos $) ' e-O 9(£ COS $) ' 

^^(^' y) = — — ' ^yi^^ y) = — — ' 

dY{x,y) dY{x,y) 

" ' ^) = 

we get the following equations of the first order in e cos 

r] = X, + Xysi + e cos + Xj,,0, ^ = {^xs + ^ cos $y,,)r/, (12) 

where 



X, = 0, Xys = Y^, = AoC (^AoL Arsech y/ZC - 1 

^ _ 2CV1 - LC ^ 2C{3LC - 2) 

Xs - ^ , Xye - — , (13) 

C (2A0L [L + AoVl - LC) Arsech - 2AoVl - LC - L{1 + Al] 



Y 

L 

Coming back to set (fTOl) and expanding the equations up to the second order in e cos $, we 
obtain the equations 

r] = Xs + Xy,i + e cos $(X, + X^.O + 17 cos^ $(X,2 + X^.a^), 

(14) 

^ = (K,, + e cos $y,e + — cos^ ^Y.^e^)ri, 
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Figure 3: Dependence of tlie y-coordinate of tlie central point yc of tlie SO on tlie perturbation 
amplitude e witli different versions of advection equations. 
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In the following, we compare numerically some properties of the SOs, generated from the 
same fixed point (xs, y^), with the main equations of motion ([8]), with linearised equations (ITOll . 
with linearized Eqs. [12] with 0(e) and Eqs. [13] with 0{e'^). Simulation shows that dependence 
of x-coordinate of the initial position of the SO on the perturbation amplitude e is practically 
the same with all the versions of the advection equations listed above. As e increases, the 
x-coordinate shifts to the west from the point Xg = 0. As to y-coordinate, it moves to the north 
with increasing e in a similar way for Eqs. [HI [ID] and [T31 but its behavior differs strongly for the 
first-order equations in e f [T^ . In Fig. [3] we demonstrate the e-dependence of the ?/-coordinate 
of the central point yc of the SO which is a central point between maximal and minimal values 
oi y. In all the = is x-coordinate of the central point. It is evident that equations 

of the first order in e cannot give a correct position of the SO. The e-dependence of deviation 
of the form of the SO from ellipticity E is shown in Fig. [1^, where E is the ratio of the area 
between the SO and an ellipse with semi- axes equal to the height and width of the SO to the 
area of the ellipse. Fig. [3)d demonstrates that the form of the SO does not change dramatically 
after linearization, but its position in the phase space cannot be found correctly in the first 
order in e. We may conclude that even with small values of the perturbation amplitude the 
effect of the second harmonic 2uj is not small. 



4 Origin and bifurcations of unstable periodic orbits 

In this section we analyze the origin of the period-4 UPOs and their bifurcations that occur with 
changing the perturbation amplitude e. Using the period-4 RM (Fig. [2^), we located all the 
period-4 UPOs and the initial positions of four trajectories for each of them. They differ by the 
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Figure 4: (a) Dependence of deviation of the SO form from ellipticity E on the perturbation 
amphtude e and (b) the form of the SO at £ = 0.2. 



type of motion of the passive particles and their length / which is a length of the corresponding 
curve between two positions on the orbit separated by the time interval 4Tq. On the cylinder 
< X < 27r, all the UPOs are closed curves whose topology may be very complicated. In 
the physical space, the UPOs can be classified as rotational ones (particles rotate in the same 
frame), ballistic ones (particles move ballistically from a frame to frame) and rotation-ballistic 
ones (particles may rotate for a while in a frame, then move ballistically through a few frames 
and change their direction of motion). 

4.1 C'^g class: western ballistic UPOs associated with the 4 : 1 west- 
ern ballistic resonance 

The shortest ones among all the period-4 UPOs are western ballistic UPOs. The particles 
belonging to the class move in a periodic way to the west along such an orbit in the 

northern separatrix layer which appears between the northern (C) and (P) regions in Fig. [T]as 
a result of the perturbation. With the help of the RM in Fig. [2^, we located two orbits 
with initial positions of 4 periodic trajectories belonging to each of them. In order to track 
out the origin of the e and h orbits, shown in Fig. [5^1 at e = 0.0785, we decrease the value of 
the perturbation amplitude e, compute the corresponding period-4 RMs, locate the C^^g orbits 
and measure their length. 

The result may be resumed as follows. The western ballistic resonance 4 : 1 with one elliptic 
and one hyperbolic points appears under a perturbation with a very small value of e. On the 
Poincare section in Fig. [Sjo at £ = 0.005, it is manifested as 4 ballistic islands along the northern 
border of the separatrix layer. Two of the islands are so thin that they are hardly visible in 
the figure. With increasing e, the size of the resonance decreases and at the critical value 
e ~ 0.016 it vanishes (see Fig. [2)3 where there are no signs of that resonance at e = 0.0785). 
The orbit, associated with the elliptic fixed point of the western resonance 4:1, loses its 
stability and bifurcates into the C^g UPO of period-4 which we denote by the symbol e. Its 
length practically does not change with increasing e (see Fig. [6]). The length / of the h orbit, 
associated with the hyperbolic point of the that resonance, changes dramatically at the point 
of bifurcation e ~ 0.0715 (see Fig. [6]) increasing fastly after this point because of appearing a 
meander and a loop on the h orbit nearby the SO. 
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Figure 5: (a) Two western ballistic period-4 UPOs at e = 0.0785 with initial positions of 
4 unstable periodic trajectories belonging to each of them. The e and h orbits were born from 
the elliptic and hyperbolic points of the northern ballistic resonance 4:1, respectively, (b) The 
Poincare section at £ = 0.005 with the northern baUistic resonance 4 : 1 consisting of 4 ballistic 
islands along the northern border of the separatrix layer. 




Figure 6: Dependence of the length / of the e and h western ballistic UPOs on the 

perturbation amplitude e. 
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Figure 7: (a) Metamorphoses of the h Cgg eastern baUistic UPO as the perturbation amphtude e 
changes: sohd hne, e = 0.0785; dashed hne, e = 0.03; and dotted hue, e = 0.01. (b) Dependence 
of the length / of the h C^g orbit on the perturbation amphtude e. 

4.2 C|g class: eastern ballistic UPOs associated with the 4 : 1 east- 
ern ballistic resonance 

The eastern balhstic UPOs of period-4 he in the southern separatrix layer between the (J) and 
(C) regions (see Fig.[T]) to the north from the jet. Particles move along such orbits to the east in 
a periodic way. The origin and bifurcations of the C^g orbits are similar to the western ballistic 
ones. They appear from elliptic and hyperbolic points of the eastern ballistic resonance 4 : 1 as 
e increases. The elliptic orbit of this resonance loses it stability and bifurcates into the UPO 
of the e-type. The hyperbolic orbit changes its topology, transforming from the bell-like curve 
at e = 0.01 to the curve with a meander and a loop nearby the SO at e = 0.0785 (see Fig. Uh). 
The length / of the h orbit in Fig. [Tb decreases with increasing e up to the bifurcation point 
e ~ 0.0275, after which it increases fastly due to the comlexification of the orbit's form (see 
Fig.CK). 



4.3 C^^ class: orbits associated with the 4 : 1 rotational resonance 

Four different orbits in the C^^ class, located with the help of the period-4 RM, are shown 
in Fig. [8]at £ = 0.0785 along with initial positions of 4 unstable periodic trajectories on each 
of them. The corresponding fluid particles rotate in the same frame along closed curves. The 
genesis of the C^^ orbits is the following. A 4 : 1 rotational resonance appears under a small 
perturbation As the amplitude e increases, its elhptic point loses stability and bifurcates 
at £ ~ 0.003 into the UPO which we denote by the symbol e. The hyperbolic orbit of the 4 : 1 
resonance h bifurcates at e ~ 0.040945 into two hyperbolic orbits, hi, h2 and an elliptic orbit 
in the center of the rotational 4 : 1 resonance. It is a pitchfork bifurcation. In Fig. [9^ we 
show by arrows a movement of initial positions of the C^^ orbits as e decreases. The pitchfork 
bifurcation point is shown as a black circle. As e increases further, the elliptic orbit loses its 
stability and becomes a hyperbolic C^^ orbit. 



4.4 C^^ class: orbits associated with the 2 : 1 rotational resonance 

The class C"^^ consists of two rotational UPOs, associated with the 2 : 1 rotational resonance, 
with elliptic and hyperbolic orbits and two periodic trajectories on each of them. The resonance 
appears under a small perturbation. At £ ~ 0.0665, the elliptic orbit of this resonance undergoes 
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Figure 8: Four 4 : 1 rotational UPOs of the C^^ class with initial positions of 4 unstable 
periodic trajectories belonging to each of them. The e and /?.i,2.3 orbits were born from the 
elliptic and hyperbolic points of the 4 : 1 rotational resonance, respectively {e = 0.0785). 




Figure 9: (a) Movement of initial positions (a;(to), y(io)) of the rotational 4 : 1 UPOs with 
decreasing the perturbation amplitude is shown by arrows. The black circle is a pitchfork 
bifurcation point at £ 0.040945 where the hyperbohc point of the resonance h bifurcates into 
three different orbits /ii,/i2 and h^. (b) Bifurcation diagram l{e). 
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Figure 10: (a) 2 : 1 rotational UPOs with initial positions of 4 unstable periodic trajectories 
{e = 0.0785). (b) Poincare section of the 2 : 1 rotational resonance at e = 0.0668 (just after 
the period-2 bifurcation). 



Figure 11: Emergency of the doubly-branched UPO, shown in Fig. [TOb. from the elliptic orbit 
of the 2 : 1 rotational resonance at the bifurcation point (black circle). Movement of fixed 
initial points on the upper /i(to) and lower h{tQ + 2Tq) branches to the bifurcation point with 
decreasing e is shown. 

a period-2 bifurcation into a period-4 elliptic orbit with 4 trajectories and a period-2 hyperbolic 
orbit with 2 trajectories. The latter is shown in Fig. [TUk at e = 0.0785. The Poincare section 
of the 2 : 1 rotational resonance at e = 0.0668 (just after the bifurcation) is shown in Fig. [TUb . 
By further increasing e, we see that the resonance 2 : 1 vanishes, and the period-4 elliptic orbit 
loses its stability and bifurcates into a period-4 UPO. To clarify this bifurcation we fix a point 
at the moment to on the upper branch of the orbit in Fig. [TOh and another point on the lower 
branch at t = to + 2To and scan their positions when decreasing the perturbation amplitude e. 
Figure [TT] demonstrates clearly that the points move to the position of the elliptic orbit of the 
2 : 1 resonance and merge with it at the bifurcation point. 

4.5 class: rotation-ballistic UPOs associated with the rota- 

tional-ballistic resonance 

The class C^g consists of four period-4 UPOs shown in Fig. W2h . We call those orbits as 
rotational-ballistic (RB) ones because the corresponding particles begin to move to the west 
in one frame, then turn to the east and travel in the southern separatrix layer to the next 
frame, fulfill one turnover in this frame and repeat their motion to the east. The genesis 
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Figure 12: Four rotational-ballistic period-4 UPOs with their initial positions in the phase space 
eh, hi, h2 and heh. The other points are the initial positions of the corresponding trajectories. 

of the C^g orbits differs from the genesis of the other classes of period-4 UPOs. Each of 
the resonances, associated with C^^, C^e, C^^ and C^'^, appears under an infinitely small 
perturbation amplitude e. Rotation-ballistic resonances of period-m and corresponding orbits 
cannot in principle appear in the flow below some critical value of e (which depends on m) 
because the width of the stochastic layers (which increases with increasing e) should be large 
enough in order that particles would have enough time to travel the corresponding distance to 
the east and west. 

The genesis and evolution of the C^g UPOs are shown schematically on the bifurcation 
diagrams in Fig. [131 In Fig. [T3k we demonstrate movement in the phase space of initial 
positions of the C^g orbits with decreasing the perturbation amplitude e. The dependence of 
the lengths of those orbits on e is shown in Fig. [T3b . The rotational-ballistic 4 : 1 resonance 
appears at a critical value of the perturbation amplitude e = ei ~ 0.040715 and manifests itself 
as 4 small islands of stability on the corresponding Poincare section. One of these islands is 
shown in Fig [T^ . If e increases further, the elliptic orbit eh in the center of the RB 4 : 1 
resonance loses its stability (it was checked by computing an accumulated error) and becomes 
a hyperbolic RB period-4 orbit. The hyperbolic orbit h of the RB resonance at e = £2 ~ 0.0433 
undergoes a pitchfork bifurcation into two hyperbolic RB period-4 orbits hi and h2 and the 
elliptic orbit heh in the center of the new period-4 RB resonance one of whose stability islands 
is shown in Fig. [TDd. Under further increasing e, the elliptic orbit heh loses its stability and 
becomes a hyperbolic RB period-4 orbit. 



5 Conclusion 

We have developed a method to track unstable periodic orbits based on computing and ana- 
lyzing local minima of the distance function d{x{tf)),y{to)). The method is not restricted to 
our simple Hamiltonian system with a periodic perturbation and is applicable to a variety of 
chaotic systems. The results, obtained with our specific model of a meandering jet, can be 
resumed as follows. 

We detected and located the UPOs of periods 1 and 4. Varying the perturbation amplitude, 
we have found bifurcations of the period-4 orbits, the UPOs of the lowest period with non-trivial 
origin. Those orbits have been grouped into five classes by their origin and bifurcations. 
(C'eb) class consists of the western (eastern) ballistic UPOs associated with the 4 : 1 western 
(eastern) ballistic resonance. C^^ (C'r^) class consists of the rotational UPOs associated with 
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Figure 13: (a) Arrows show movement of the initial positions of the C^g orbits with decreasing 
the perturbation amphtude e. At the point e = Si, there appears a RB 4 : 1 resonance with 
the elhptic orbit eh which loses its stability with increasing e. The hyperbolic orbit of the RB 
resonance bifurcates aX e = 62 > £1 into two hyperbolic period-4 UPOs hi and h2 and the 
elliptic orbit heh which loses its stability with further increasing e. (b) Bifurcation diagram 
l{e). 
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Figure 14: (a) One of the stability islands of the RB 4 : 1 resonance appearing after the first 
bifurcation at £ = 0.04072 > Si with the elhptic point eh in its center, (b) One of the stabihty 
islands of the RB 4 : 1 resonance appearing after a pitchfork bifurcation at £ = 0.0435 > £2 
with the elliptic point heh in its center. 
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the 4 : 1 (2 : 1) rotational resonance. C^g class consists of specific rotational-ballistic UPOs 
associated with 4 : 1 rotational-ballistic resonance. We would like to stress the following. Ro- 
tational and ballistic resonant islands are well known objects in the phase space of Hamiltonian 
systems [2Tj. As far as we know, rotation-ballistic resonant islands of stability have not been 
found before. They are expected to appear as well in other jet-flow models (kinematic and dy- 
namic ones) under specific values of control parameters, and it would be interesting to search 
for them in the models introduced in Refs. [a[3lll[5l[6l[ia[l3l[T3[l5l[l6l[I7l[l8l[l9l[2^ It 
is expected that boundaries of the rotational-ballistic islands are specific dynamical traps that 
should affect transport and statistical properties of passive particles and may be treated by the 
methods developed in 

In addition, by linearizing the advection equations, we have studied the properties of the 
period- 1 saddle orbits playing a crucial role in chaotic transport and mixing of passive particles. 
The UPOs with larger values of the period can be classified into three big groups: ballistic, 
rotational and rotation-ballistic ones. The origin and bifurcations of the orbits in each class 
can be studied in a similar way, but it may require larger computational efforts. 

The results obtained do not depend critically on the exact form of the streamfunction and 
chosen values of the control parameters. The UPOs in other kinematic and dynamical models 
of geophysical jets, known in the literature [121 IISI [11 [ISl |16l IIIl [TBI [ISl EO] , can be detected, 
located and classified in a similar way. The UPOs form the skeleton for chaotic advective mixing 
and transport in fluid flows, and the knowledge of them (at least, lower-period ones) allows us 
to analyze complex, albeit rather regular, stretching and folding structures in the flows. 
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